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Based on a supersymmetric Yukawaon model with 0(3) family symmetry, possible forms of quark 
and lepton mixing matrices are systematically investigated under a condition that the up-quark mass 
matrix form leads to the observed nearly tribimaximal mixing in the lepton sector. Although the 
previous model could not provide a good fitting of the observed quark mixing, the present model can 
give a reasonably good fitting not for lepton mixing but also the quark mixing by using a different 
origin of the CP violation from one in the previous model. 
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t^H . I. INTRODUCTION 

One of the current subjects of the particle flavor physics is to understand quark and lepton masses and mixings. 
The investigation of them, even if it is phenomenological, will provide a promising clue to new physics. The so-called 
• "tribimaximal" mixing observed in the neutrino mixing [l], Q is very suggestive of a fundamental law in the lepton 
Oh| sector. Usually, the observed tribimaximal mixing has been explained by assuming a discrete symmetry Q. 

Meanwhile, as a neutrino mass matrix model without assuming any discrete symmetry, an unfamiliar model 0, HI 
has been proposed by using a seesaw type neutrino mass matrix M v = m^M^m^. In this model, the Dirac neutrino 
. mass matrix mrj is given by mr> oc M e (M e is a charged lepton mass matrix) and the Majorana mass matrix Mr of 
the right-handed neutrinos is given by 

*> '■ M R oc Ml' 2 M e + M e M^ 2 . (1.1) 

<n : 

t-H , Here Mu" 1 is defined by M^ 2 (M^ 2 ) T = M u (M u is an up-quark mass matrix with a symmetric form Mj = M u ). 
The model d can lead to a nearly tribimaximal mixing by assuming suitable up-quark mass matrix as we give a short 
review in the next section. The model has only four parameters: one (£„) is in the neutrino sector, and one (a u ) is 

t-H ■ in the up-quark sector, and two {a.de lad ) in the down-quark sector. (Here, we consider that the charged lepton mass 
, (m e , to m , m T ) are known parameters, and we do not count these parameters as adjustable parameters.) This model 
leads to excellent fitting for up-quark mass ratios m u /m c and m c /mt and neutrino mixing parameters sin 2 28 a t m , 
tan 6 so i ar and |J7 e 3| 2 , only by adjusting the two parameters a u and £„. On the other hand, for down-quark sector, 
the fitting is not so excellent, especially, the predicted values of \V u b\ and \Vtd\ are somewhat large compared with the 

k>( | observed values as far as we use parameter values which can give reasonable values for the observed down-quark mass 
. ratios. 

The purpose of the present paper is to investigate an improved version of the above model and to search sys- 
tematically for parameter values which can give reasonable quark mass ratios, quark mixing parameters (Cabibbo- 
Kobayashi-Maskawa (CKM) mixing matrix) and neutrino mixing parameters. In Sec. II, we will show that the four 
parameter model cannot have reasonable parameter region consistent with four quark mass ratios, three neutrino 
mixing parameters, and four CKM mixing parameters. In Sec. Ill, we propose a revised model and give parameter 
fitting for 11 observables. (In the present model, we do not discuss the observed value R = Am^ olar / Am atm for 
neutrino masses, because we can always have an additional one parameter which inevitably appears in the model and 
affects only the mass ratios R, but does not affect neutrino mixing and observables in the quark sector.) Finally, 
Sec. IV is devoted to the summary and discussions. 
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II. SUPERSYMMETRIC YUKAWAON MODEL 

In this section, we give a short review of a quark and lepton mass matrix model Q based on the supersymmetric 
Yukawaon model, because, in this paper, we propose a revised version of this model. 

In the Yukawaon model, we put the following assumption: 
(i) We consider that the Yukawa coupling constants are effectively given by 

yf = yf^ (2-1) 

where (Yf) (f = u, d, e, v and so on) are vacuum expectation values (VEVs) of new scalars Yf with 3x3 components 
of 0(3) family symmetry and A is an energy scale of the effective theory. (For the time being, we assume A ~ 10 14 ~ 15 
GeV.) Therefore, the would-be Yukawa interactions are given by 

Hy = E ^(Y a ) ijqj H a + y fdi{Y d ) l3q] H d 



+ J2 ^i{Y»)irtH u + ^(Y^Hd + H.c. + VrvSOTh)^, (2.2) 
where q and £ are SU(2) l doublet fields, and f c (f — u, d, e, v) are SU(2) l singlet fields. 

(ii) In order to distinguish each Yf from others, we assume a \J(l)x symmetry (i.e. "sector charge") in addition to 
the 0(3) symmetry, and we have assigned U(l)x charges as Qx(Yf) = Xf, Qx(f c ) — —Xf and Qx(Yr) = 2x v . (The 
SU(2)i doublet fields q, £, H u and Hd are assigned to sector charges Qx — 0.) 

(iii) For the neutrino sector, we assume Qx{v c ) — Qx(e c ), so that the Yukawaon Y e can also couple to the neutrino 
sector as (tY e v c )H u instead of {£Y u v c )H u in Eq.(2.2). Therefore, we can change the above model into a model without 
Y u . Hereafter, we read Y v in Eq.(2.2) as Y e . Besides, we can have a term 



E 



Vfv${Y e ) ik {Y e ) kj v<;, (2.3) 



in addition to the right-hand side of Eq.(2.2), because Y e Y e has the same U(I)x charge as Yr, i.e. Qx = 2%e- 
Although this term (2.3) leads to an additional neutrino mass term, the term does not affect neutrino mixing as 
far as the neutrino mass matrix M v is real, 1 because of M v oc Y e [(- • • ) + Y e Y e ] _1 Y e = [Y e _1 (- • • )Y ( T 1 + 
(iv) We give a superpotcntial W which is invariant under 0(3) family symmetry and XJ(l)x symmetry, and solve 
supersymmetric (SUSY) vacuum conditions. As a result, we obtain VEV relations among Yukawaons. 
For example, we have assumed the following superpotential 

W e = A e Tr[$ e $ e O e ] + MeTr[Y e 6 e ] + W*. (2.4) 

Here we have assumed Qx($ e ) = \Qx{Y e ) — — \Qx{®e) and the term W$ has been introduced in order to determine 
a VEV spectrum ($ e ) completely. Then, from a SUSY vacuum condition 

dW 

— = Xe^e + VeYe = 0, (2.5) 
OBe 

we obtain a VEV relation 

(Ye) = --<$e><$e>. (2.6) 
Me 

The VEV value ($ e ) is derived from the term W$ (for example, see Refs.0-@)- However, for simplicity, in this paper, 
we use the observed values of the charged lepton masses straightforwardly for the VEV value as given by 

($ e ) e = diag(vi,V2, V3) oc diag(y / ml, V"V> V™t)- (2-7) 



When R T MiR = Di = diagonal, the inverse matrix M-t 1 is also diagonalized as R T Mj 1 R = D 1 1 by the same orthogonal trans- 
formation matrix R; When we take M = Mi + mol, M is diagonalized as R T MR = Di + mol, so that we can diagonalize M~ 1 as 
R T M~ 1 R = (R T MR)- 1 = (D + mol)" 1 . 
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In other words, we have assumed the ad hoc relation (2.7), the derivation of which is not discussed in the present 
paper. Hereafter, for counting a number of "adjustable" parameters, we do not include Vi in the number. Here, the 
notation (A) f denotes a form of a VEV matrix (A) in the diagonal basis of (Yt) (we refer to it as / basis). The 
scalar O e does not have a VEV, i.e. (O e ) = 0. Therefore, terms which include more than two of e do not play any 
physical role, so that we do not consider such terms in the present effective theory. [Hereafter, we will denote fields 
whose VEV values are zero as notations 0a (A — u, d, ■ ■ ■ ).] 

Next, for the purpose of the comparison of our new model with the previous one, we give a short review of quark 
and lepton mass matrix forms of the previous model discussed in Ref. [5| . The explicit form of the superpotential for 
the previous model is given in Ref.@. That, for the new model, shall be given in the next section. 

In the previous model [5], the quark mass matrices, i.e. (Y u ) and (Yd), are given as 

M„oc(y u )oc <$„)<$„), (2.8) 
($ u ) e cx ($ e ) e ({E) e + a u (X) e ) ($ e ) e , (2.9) 

M d cx (V d ) e cx ($ e ) e ((^) e +a d e Md (V) e )($ e ) e , (2.10) 

respectively. (For convenience, we have changed the definitions of a u and ad from those in Ref. H.) Here, (X) e and 
(E) e are given by 

x (1 1 1\ /I 0\ 

(X) e = Xv x = -\ 1 1 1 \v x , (E) e = lv E = 1 \v E . (2.11) 
6 \ 1 1 1 / \ 1 / 

(Here, the VEV form (X) e breaks the 0(3) flavor symmetry into S3.) Therefore, we obtain quark mass matrices 

M*/ 2 oc Mg/ 2 (1 + a u X) M*/ 2 , M d oc M X J 2 (l + a d e la "X) M 1 / 2 , (2.12) 

on the e basis. Here, we have redefined the coefficients a u vx/vE and a d vx/vE in Eqs.(2.9) and (2.10) as a u and a d , 
respectively. Hereafter, for numerical estimate of a u an d au , we use the definition of those in Eq.(2.12). (This quark 
mass matrix form (2.12) has first been proposed in Ref. }lQj| as a "democratic universal seesaw mass matrix model".) 
Note that we have assumed that the 0(3) relations are valid only on the e and u bases, so that (Y e ) and (Y u ) must be 
real. [The VEV matrix (3> u ) must satisfy the relation (2.9) on the e basis, while must also satisfy the relation 
(Y u ) oc ($u)u(®u)u on the u basis. However, for the down-quark sector, such a condition is not required, because (Yd) 
is given by Eq.(2.10) only on the e basis.] 

A case a u ~ —1.79 can give a reasonable up-quark mass ratios \Jm u \jm u i — 0.043 and ^/m U 2/m U 3 — 0.057, which 
are in favor of the observed values fll|j 




772 

= 0.045^-^1 ./—= 0.060 ±0.005, (2.13) 
m c ' V m t 

at /j, = mz- 

In this paper, we will carry out parameter- fitting at fx = mz, because we interest in the mixing values at fj, = mz- 
Exactly speaking, fitting for the mass ratios must be done at fi = A ~ 10 14 GeV. However, at present, our model does 
not intend to give so precise predictions of the quark mass ratios. For example, we know [ll[ yj m u /m c = 0.046lo oi2 
and yjm c jra t = 0.0511q'qo6 even at /i = 2 x 10 16 GeV (tan/3 = 10). Even in \fm c /mu the discrepancy is smaller 
than 20%. Besides, the mass values are dependent on the value of tan/? in the SUSY model. Therefore, for simplicity, 
in this paper, we will carry out the parameter-fitting at /1 = mz- 

On the other hand, in the neutrino mass matrix 

M v oc (Y e ) e (Y R )^ 1 (Y e ) e , (2.14) 
the Majorana mass matrix of the right-handed neutrinos (Yp) e is given by 

(Y R ) CX (Y e ) e (P u ) e ($ u ) e + ($ u ) e (Pu)e(Y e )e + Zv((Pu)e(Y e )e($u)e + (®u)e (Ye) e (P„) e ). (2.15) 

Here, we have introduced a new field P u with a VEV 

(P u ) u ocdiag(+l,-l,+l), (2.16) 
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FIG. 1: \Vij\f\Vij\obs versus a phase parameter a d . (a) |Ka|/0.2252 and | V c6 |/0.0406; (b) | V ub |/0. 00389 and |V td [/0.0084. The 
parameter a u in the up-quark sector is fixed at a u = —1.79. The five curves represent a d = —1.30, —1.40, —1.50, —1.60 and 
— 1.65, respectively. 



in order to make "effective" eigenvalues of (<&„)„ positive, because the eigenvalues of (<&„)„ = v u %, v us) give signs 
(+, — , +) for the parameter value a u ~ —1.8. (The field P u has been introduced from a phenomenological reason. If 
the factor (2.16) were absence [i.e. Yr were given by (Yr) oc {Y e ) e (<fr u } e + ($„) e (Y e ) e ], we could not give the observed 
maximal mixing sin 2 29 atm — 1 [H for any values of the parameters a u and The reason for the £„ term in Eq.(2.15) 
is as follows: When we consider a term Y e P u Q u + $ u P u Y e we must also consider an existence of P u Y e $> u + & u Y e P u 
H, because they have the same U(l)x charges. The results at a u ~ —1.78 are excellently in favor of the observed 
neutrino oscillation parameters sin 2 29 atm = 1.00_o 13 [H and tan 2 8 so i ar — 0.469+o'o4i by taking a small value of 
16/1, tu = +0.005 or f„ = -0,0012. 

Thus, the model in Rcf.[5] can successfully fit two up-quark mass ratios and three neutrino mixing parameters 
only by the two parameters a u and £„. On the other hand, the fitting of six observable quantities (two down-quark 
mass ratios and four CKM mixing parameters) only by two parameters ad and ad given in Ref. [f| are not in excellent 
agreement with the observed values. Especially, the predicted values of \V u b\ and \ Vtd\ are considerably larger than the 
observed ones. We find from a systematical parameter search that this is not due to incompleteness of the parameter 
search, but plausible values of the CKM mixing parameters cannot be obtained even if we abandon the fitting of the 
down-quark mass ratios. 

Considering the success in the up-quark and neutrino sectors, we do not change the model for the up-quark and 
neutrino sectors. We fix the parameter values as a u ~ —1.8. The observed values of the down-quark masses are as 
follows [13 

— = 0.05312-jgL — - 0.019 ± 0.006, (2.17) 

at (J, = mz- We consider that the mass ratio md/m s may be sensitive due to an unknown effect of a minor change of 
the model, so that, for the time being, we disregard the fitting of md/m s and concentrate on the fitting of m s /mb- 
Although a parameter value ad ~ —16 can give a reasonable prediction of m s /mb, the solution cannot give reasonable 
predictions of the CKM mixing parameters, we rule out this solution ad ~ —16. We find that there are another 
solutions of ad in a range —ad — 1.3 — 1.7, which can roughly give m s /mb = 0.1 — 0.4. The solutions have a 
possibility that they can give reasonable values of the CKM mixing parameters. Therefore, we investigate the case 
with —ad = 1.3 — 1.7 in detail. 

The results are shown in Fig.l, where the predicted values \Vij\ versus the phase parameter a.d are given in the unit 
of the observed values [l2j |V^j| 6 S 



\V us \ obs = 0.2252 ± 0.0009, |V cfc | o6s = 0.0406 ± 0.0013, 
\Vub\obs = 0.00389 ± 0.00044, \V td \obs = 0.0084 ± 0.0006. 



(2.18) 
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Here, we have illustrated the behaviors of \Vij | for the range ay = 0° — 180°, because the behaviors for ad = 360° — 180° 
are just the same as that for ad = 0° — 180°. As seen in Fig. 1(a), in order to obtain a reasonable value of \V c b\, we 
must choose a value of ad smaller than a<j ~ 10°, and also a value of a<j smaller than — ad ~ 1.5. However, from 
Fig. 1(b), we can conclude there is no solution for a reasonable value of \V US \ for any values of ad and ad even at the 
cost of the fitting of down-quark mass ratios. Therefore, in the next section, we proposed a revised model for quark 
mass matrices keeping the model for the neutrino sectors. 



III. PHENOMENOLOGY OF QUARK MASS MATRICES 



We present an explicit form of the quark mass matrices in our new model. In this paper, we put the following 

1 /2 

assumptions for a phcnomcnological forms of quark mass matrices M u [M u ) and Md'. 

(i) Differently from the previous model |5J, we regard that not only (Y u ) (also but also (Yd) are real, i.e. 

ad — in Eq.(2.10). Instead, we consider that CP violation in the quark sector originates in a phase matrix 
(Pd)e = vpdPd = wprfdiag(e 1 * 1 , e l< ^ 2 , e 1 ^ 3 ) which does not affect the down-quark mass ratios, but does only the CKM 

1/2 

mixing. Namely the quark mass matrices M u and Md is given by 

M X J 2 oc M e 1/2 (1 + a u X)Ml /2 + (M e 1/2 A/ e 1/2 (l + a u X) + (1 + a„A)M e 1/2 M e 1/2 ) + m 0u l, (3.1) 



M d cx Pd 



M X J 2 (\ + a d X)Ml' 2 + & {m\I 2 M\I 2 (\ + a d X) + (1 + a d X)M^ 2 M^ + m od l] P d , (3.2) 



so that the CKM matrix V is given by V — R^PdRd, where R u and Rd are defined by R^Mu R u = 
diag(+ % /m^, -y/m~ c , +- v /m^) (for a u ~ -1.8) and Rj(PjM d P^)R d = diag(m d , m s , m b ), respectively. 

(ii) Similarly to Eq.(2.15), we assume the £ g terms which originate in the reordering of the fields with the same 
quantum numbers. 

(iii) Since only two of the three phase parameters fa, fa and fa in the phase matrix P d = diag(e i( ^ 1 , e** 2 , e 4 * 3 ) are 
physically independent parameters. For convenience, we take fa — 0. 

(iv) It is better that the parameter number is as few as possible. We consider that the first term is dominant in 
Eq.(3.1) [and also Eq.(3.2)], and we will consider £ g and mo, terms as the need arises. As seen later, we can do fitting 
without £d and mou terms. 

Of course, we consider that these relations are derived from SUSY vacuum conditions for a given superpotential 
W . However, prior to investigating the superpotential form, from a phenomenological point of view, we would like 
to investigate whether there is a possible parameter region or not in the present model. A Yukawaon model for the 
phenomenological forms (3.1) and (3.2) will be discussed in the next section. 

1/2 

Since the mass spectra M u (a u ) and Md(a,d) have the same behavior for the parameter a q (a u and ad), we illustrate 
the mass spectra versus a q in the limit of ^ q — in Fig. 2. (The mass values in Fig. 2 read (\/|m u |, ^/|m c |, x/J^tj) and 
(|m<j|, |m s |, \rrib\) for the up- and down-quark sectors, respectively. 

In the present model, too, the model for the up-quark sector and neutrino sector is essentially unchanged from 
the previous model [5j except for the £ M term given in (3.1). For reference, in Fig. 3, we illustrate the up-quark mass 
ratios ^Jm u /m c and yj m c /mt versus a u and £„. As seen in Fig. 3, there are two set of the solution (a«,£«) [regions 
(i) and (ii) illustrated in Fig. 3] which can give reasonable up-quark mass ratios. However, the region (ii) cannot give 
reasonable CKM mixing parameters. Hereafter, by taking fitting of neutrino mixing parameters into consideration, 
too, we will take a u = —1.764 and £ u = 0.0070 in the region (i). The choice of (a u ,£ u ) — (—1.764,0.0070) can give 
up-quark mass ratios 

& = 0.0619, ,/^ = 0.0559. (3.3) 

V m c V m t 

In model building of the down-quark sector, we give the down-quark mass ratio m s /mi preference rather than 
TOd/TO s , because it is not so difficult to adjust the ratio md/rn s without affecting the CKM parameter fitting as we 
demonstrate later. In Fig. 4, we illustrate behavior of m e /mb versus ad- As seen in Fig. 4, there are three regions 
which can give reasonable mass ratio m s /mb. However, the regions (ii) and (iii) cannot give reasonable CKM mixing 
parameters. (The region (ii) corresponds to a parameter region adopted in the old model Q.) Hereafter, we will show 
the region (i) (i.e. ad ~ —17) in detail. 

Next, we investigate possible parameter regions which can give reasonable CKM mixing parameters. We take \V US \, 
\V u b\, \V c b\ and |Vtd| as four independent parameters in the CKM matrix. In Fig. 5, we illustrate allowed regions in 
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FIG. 2: Eigenvalues |mi |, \mi\-, [ma| versus a parameter a in a mass matrix M = JVfi' /2 (l+ae lc *X)Mi' /2 with a = (thick curves). 
For reference, a case with a = 20° is also illustrated by thin curves in the figure. (|fni|, |wia|, \rr13 1) read (y / |m 1 i||, ^/|m c |, -\/|mt|) 
and (|md|, |m s |, |m(,[) for up- and down-quark sectors, respectively. Numerical values of the eigenvalues are given in a unit of 
(m e + m M + m T ). 
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FIG. 3: Allowed region in the a u - £ u plane obtained from the up-quark mass ratios \J m u /m c and \/ m c /mt- The shaded areas 
are consistent with the observed values given in Eq.(2.13). 

the 4>i ~ 4>2 plane obtained from \Vij\ with |14,| b s under a u = —1.764, = 0.0070 and = —16.6 whose values 
are obtained from global best fit. As seen in Fig. 5, the value of (f>2 — 180° is in favor of the observed CKM mixing 
parameters. The case with (j>2 — 180° is also illustrated in Fig. 6. It is interesting that \Vtj\ take their minimum at 
4>i ~ 180°. From Fig. 5 We find that (/)% ~ ±16° + 180° is in favor of the observed CKM mixing parameters. 
In conclusion, our best hit parameters are 

a u = -1.764, Cu = 0.0070, a d = -16.6, fa = 196.0° (164.0°), fa = 181.5° (178.5°) (3.4) 

together with ^ = 0, and then we obtain the predicted CKM mixing parameters 

|K» | = 0.2259, \Vd,\ = 0.04141, \V ub \ = 0.00418, \V td \ = 0.00854. (3.5) 

However, the parameter value a c i = —16.6 gives considerably small value of mj/m s , i.e. md/m s = 0.00358. In order 
to correct this wrong value, we must take mod with a non-zero value. By taking a value 



mod/mo — —0.0061, 



(3.6) 



7 



ms 
mb 



0.14 






- 


0.12 






; 


0.10 








OS 








0.06 


; ( 




(ii) ( 


ii) 


0.04 










0.02 










0.00 







-40 -20 20 40 

FIG. 4: ad dependence of m a /mb- The dotted lines show the observed down-quark mass ratio m s /nid = 0.019 + 0.006 and 
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FIG. 5: Allowed regions in the <j>i - <j>2 plane obtained from the CKM mixing parameters \Vij\. Shaded areas are consistent 
with the observed values |Vij|obs in Eq.(2.18) . The parameter values of a u , £u and ad are chosen as a u — —1.764, = 0.0070 
and ad = —16.6 respectively. The star (*) indicates the best fit points [see Eq.(3.4)]. 



where tuq = m e + m M + m T , we obtain the reasonable down-quark mass ratios 

= 0.0529, — = 0.0231, (3.7) 

m s m b 

without affecting the CKM mixing parameters. 

On the other hand, for neutrino mixing parameters, the model is essentially the same as before. In Fig. 7, we 
illustrate £„ dependence of the neutrino mixing parameters. As seen in Fig. 7, the model can predict sin 2 29 atm — 1 
and 1 1/ 13 1 2 ~ independently of £ v . The value of £„ is determined from the observed value tan 2 9 so iar = -4:57^^1. 
The value of £„ and the neutrino mixing parameters are listed in Table I. 

Let us summarize above phenomenological considerations for the mass matrices for quarks and neutrinos. By taking 
the phenomenological considerations £ f j = and mrju = into consideration, we have adopted the quark mass matrices 

1/2 

M u ' and Md given by 



Ml' 2 cx M e V2(i + a u X)Ml' 2 + (Af e 1 / 2 M e 1 / 2 (l + a u X) + (1 + a u X)M^ 2 M^ , 



(3.8) 
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FIG. 6: <f>i dependence of the CKM mixing parameters \Vij\ in the case (f>2 = 180°. In the case of <f>2 = 180°, 4>i = 191° and 
169° are in favor of all the observed | Vij | ea; p (2cr) . 




FIG. 7: dependence of the neutrino mixing parameters tan 2 6 so i ar (solid curve), sin 2 26 a tm (dashed curve) and | tT" 13 1 2 (dotted 
curve). The up-quark mass matrix parameters are chosen as a u = —1.764 and £ u = 0.0070. 



M d cx P d 



M^ 2 (l + a d X)M^ 2 + m od l 



(3.9) 



On the other hand, the neutrino mass matrix is given by Eqs.(2.14)-(2.16). By using these mass matrices with the 7 
free parameters, a u , £ u , a d , <fii, fa, mod, and we have searched systematically for the parameter values which can 
give reasonable 4 quark mass ratios, 4 CKM quark mixing parameters, and 3 neutrino mixing parameters. (Although 
the values of (m e , m^, m T ) play an essential role in the present model, we have fixed those to the running mass values 
at ji = mz, so that we do not count those as free parameters.) 
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tan Qsolar 


sin 2 20atm 


\Ui 3 \ 2 


+0.00031 


0.457 


0.999 


2.56 x 10" 4 


-0.00102 


0.457 


0.998 


2.74 x 10 -4 



TABLE I: Input value of and predicted values of the neutrino mixing parameters. The up-quark mass matrix parameters 
are chosen as a u — —1.764 and £ u = 0.0070 which can give reasonable up-quark mass ratios. 

IV. SUPERPOTENTIAL 

In this section, by taking the phenomenological results with £ d = and mo« = in the previous section into 
consideration, we discuss a possible form of the superpotential W assuming an 0(3) family symmetry. Since we 
consider the effective theory with A ~ 10 14 GeV, at present, it is not our chief concern whether 0(3) is local or global. 
For the moment, we assume that 0(3) is global. It should be noted that the massless states are harmless because A 
takes an extreme large value A ~ 10 14 GeV (l3j . Under the 0(3) family symmetry and conservations of \J(l)x and 
R charges given in Table II, we obtain the following form of W: 

W = W e + W R + W u + W d , (4.1) 
W e = M e Tr[y e e e ] + A e Tr[$ e $ e O e ], (4.2) 



Wr = ^{Y R Q R } + {Tr[(Y e P u $ u + $ u P u Y e )® R ] + ^Tr[(P„y e $ u + $„y e P„)0 fi ]} , (4.3) 

W u = M„Tr[y u 6„] + A u Tr[$ u $ u O u ] 
+ ^ {Tr[$ e S u $ e Q' u } + £„Tr[(3> e $ e S u + S u $ e $ e )G' u }} , (4.4) 

W d = ^Tr[P d Y d P d G d ] + ^Tr[$ e S d $ e O d ] + n od Tr[E od Q d ], (4.5) 

where, for convenience, we have denoted linear combinations of fields E q and X q as S q = E q + a q X q (q = u, d) and 
(E q ) = v Eq l and (X q ) = v Xq X. 

Among the SUSY vacuums which are derived from the superpotential (4.1), we take only a vacuum with (0 e ) = 
(Br) = (Q u ) = (Q' u ) = (e d ) = 0. Therefore, we can obtain VEV relations (2.6), (2.15), (3.1) and (3.2) from SUSY 
vacuum conditions 8W/dO e = 0, dW/d&R — 0, dW/d<d' u = and dW/dO d — 0, respectively. Since other conditions, 
for example, dW/dY e = 0, and so on, inevitably contain a field &a (A = u,d, ■ ■ ■ ), they cannot play effective roles 
in the VEV relations. (Although we did not give an explicit form of ^(^e) in the present paper, we assume that 
W($ e ) also contains O fields. For the form of ($ e ), Eq.(2.7), we will use the charged lepton mass values at /i = mz- 
) One of merits to introduce such O fields is that we do not need to consider contributions from higher dimensional 
terms with the form (Tr[- ■ • Oa])™ (n > 2), because dW/d&A from such a higher dimensional term always contains 
Ga more than one, so that such a term becomes vanishing. 

Let us emphasize a role of the R charges: By assuming the R charge conservation with the R charge assignment 
given in Table II, we can forbid all of higher dimensional terms with (1/A)" (n > 2) except for the terms given 
by Eqs.(4.2) - (4.5). (However, for this purpose, we must assume that our Kahler potential is given by a canonical 
(minimal) form.) We also note that if we assume U(l)x only, the assignments can allow unwelcome terms in the 
superpotential, for example, Tr[S u P u ], Tr[$ u P u O e ], and so on. Such terms can be forbidden by assuming suitable R 
charge assignments. For example, when we take P-charges as R(£) — 1— r, R{e c ) — R{v c ) — R(q) — R(u c ) — R(d c ) = 1 
and R(H U ) = R(H d ) = 0, we can forbid the terms Tr[S u P u ] and Tr[$ u P u e ] by taking R charges of other fields as 
given in Table II. 

In the phenomenological study in Sec. Ill, the VEV values of 4> e play an essential role in evaluating the predicted 
values. Although it has been tried to build a model @,|^| which gives VEV spectrum (2.7), it is not clear whether such 
a model can be applicable or not to the present model straightforwardly. In this paper, we do not give a superpotential 
form which can lead to the VEV spectrum (2.7). We have just assumed the VEV value given by Eq.(2.7), where we 
have used the values of charged lepton masses at the scale [i = mz- 
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Fields Y e $ e 9 e Y R Or Y u Pu 



Eu,X u Yd Pd Qd Ed,Xd Ead 



Qx «Ce 2"^ e Xu ~^Xu Xe ^Xu 

R charge r |r 2 - r 2 



2 



\x u TjX u — X e Xd Xp —(Xd + 2xp) Xd + 2x p — X e Xd + 2xp 

2 -r 2 -r r 



TABLE II: U(l)x charges of the Yukawaons. For the time being, we assign different charges for the fields E u and Ed (X u and 
Xd) by assuming that those are different fields. For R charges, see text. 

Also, so far, we have not given superpotential forms which lead to VEV matrices (E q ), (X q ), (P u ) and (Pd)- In 
general, any Hermitian VEV matrix (A) can be obtained from a superpotential 



[However, we must assign a U(l)x charge — 3Qx(A) for the coefficients Xi (i = 1,2,3).] For example, (Ai, A 2 ,A 3 ) = 
(1/6, —1/2, 1/3) (i.e. W = detA) gives the form (X) given in Eq.(2.11). However, this method is not applicable to 
the form (Pd), because the VEV matrix is not Hermitian. In this paper, we have assumed these ad hoc VEV forms. 



In conclusion, we have proposed a phenomenological quark and lepton mass matrices based on a Yukawaon model. 
In Sec. II, we have demonstrated that the previous model in which the CP violation originates only in the complex 
parameter ad, cannot give reasonable CKM mixing values even if at the cost of the quark mass ratios. Differently 
from the previous model, in the present model, CP violating phases are introduced in the phase matrix Pd given 
in Eq.(3.2). In the up-quark sector, we have considered the £ u term in Eq.(3.1). This comes from the fact that the 
terms with another order of the fields, $ e $ e 5 u + S u & e $ e , cannot be, in general, forbidden compared with the order 
of $ e 5 u $ e because of the same U(l)x charges. A similar situation have been assumed in the neutrino sector, too, 
i.e. the £„ term in Eq.(4.3). (The values of and £ v are very small.) In contrast to those sectors, in the down-quark 
sector, we have not considered such a term as well as an additional term PdPdYd + YdPdPd corresponding to PdXdPd 
in Eq.(4.5). This is a result from the phenomenological study, and the theoretical reason for the absence is unknown 
at present. Also we note that the phenomenological fit requires the mod term added to Eq.(3.2), but it does not need 
an rriQu in Eq.(3.1). 

Our numerical conclusions from the present systematical study is summarized in Figs. 2-7. Especially, as seen in 
Fig. 7, the results sin 2 29 atm — 1 and | CT^ X3 1 2 < 0.005 are insensitive to the value of the parameter £„. In other words, 
if Wiz\ 2 ~ 0.01 (the possibility was pointed out by Fogli, et al. [14j ) is established experimentally, the present model 
will be ruled out, or it will need a drastic revision. 

We have been able to obtain reasonable parameter fitting not only for the observed lepton mixing but also for the 
observed quark mixing. However, the model still includes ad hoc assumptions. We consider that it is important to 
clarify what parts are problems to get a good fitting of the data for the next step of the investigation. Our model 
building will proceed step by step. 
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